Transport model study of the m^-scaling for A, K, and 7r HBT-correlations 
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Based on the microscopic transport model UrQMD in which hadronic and string degrees of free- 
dom are employed, the HBT parameters in the longitudinal co-moving system are investigated for 
charged pion and kaon, and A sources in heavy ion collisions (HICs) at SPS and RHIC energies. In 
the Cascade mode, Ro and the Rl at high SPS and RHIC energies do not follow the wr-scaling, 
however, after considering a soft equation of state with momentum dependence (SM-EoS) for formed 
baryons and a density-dependent Skyrme-like potential for "pre-formed" particles, the HBT radii of 
pions and kaons and even those of As with large transverse momenta follow the m-r-scaling function 
R = 3/ y/rnr fairly well. 

PACS numbers: 25.75.Gz,25.75.Dw,24.10.Lx 
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In order to explore the properties of the fireball which 
is created just after the collision of two energetic nu- 
clei, the Hanbury-Brown-Twiss interferometry (HBT) 
technique has been used and developed in astro- and 
nuclear-particle physics for about half a century. It is 
well-known that the HBT technique can provide impor- 
tant information about the spatio-temporal structure of 
the particle emission source (the region of homogene- 
ity). In [![, a non-trivial transition in the excitation 
function of spatio-temporal characteristics of the source 
was proposed. Meanwhile, the AGS, SPS and RHIC ex- 
periments (with nucleon-nucleon center-of-mass energies 
^Jsmn from about 2.5 GeV up to 200 GeV for heavy ion 
collisions) have been stimulating the HBT related inves- 
tigations further into another golden era. Although the 
experiments have discovered quite a few exciting hints 
for a new phase of matter - maybe a quark gluon plasma 
(QGP) (see, e.g., Refs. [3,3), a plethora of puzzling phe- 
nomena also came out. Examples related to HBT are, 

1. the discovered that there is no obvious peak/ valley 
of the HBT quantities over the whole beam energies 
from SIS, AGS, SPS, up to RHIC as suggested in 

a. 

2. that hydro-dynamic as well as transport (cascade 
mode) calculations show that the calculated ratio of 
the HBT radii in outward and sideward directions is 
higher than the experimental data, which is known 
as the HBT time-puzzle ("t-puzzle") [I S, i, 0- 
Forpossible solutions of the "t-puzzle" we refer to 

Hi on mill Hi- 
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3. mainly due to the resonance decay contribution of 
the particle production as well as the effects of 
strong and electromagnetic forces, the correlation 
function of identical particles deviates from a Gaus- 
sian type 

A "direct" study of the so-called "non-Gaussian ef- 
fect" — the emission source image techni que- — has been 
quickly improving in recent years fl5l []^.1l7l [l8j. How- 
ever, in order to theoretically investigate previous data 
and due to the fact that it gives a leading-order ap- 
proximation to the real shape of the homogeneity re- 
gion, the Gaussian parameterization is still important. 
During the process of Gaussian parameterization, it is 
known that on both experimental and theoretical sides, 
a proper potential modification of the final state interac- 
tion (FSI) after freeze-out should be analyzed before any 
work is done because of the intrinsic physical character- 
istic of correlated pairs. For example, charged pion-pion 
interferometry is only relatively weakly affected by the 
Coulomb and nuclear potentials due to its small mass 
and collision cross section, while for charged kaon-kaon 
case, the Coulomb modification should be taken into ac- 
count in analysis [3, Ell • For baryon-baryon correlation, 
a proper nuclear potential modification has to be consid- 
ered. In the absence of FSI and from a hydro-dynamic 
point of view in which the freeze-out of particles is flow- 
dominated, the HBT radii are predicted to decrease with 
l/^/rriT (mr-scaling), txit is the transverse mass of the 
observed particle-pair, independent of the particle species 
[20l ]. This prediction has been probed by recent exper- 
iments with energies from AGS to RHIC. In Ref. [21j . 
a ray-scaling expression has been suggested at 

Eiab = 158 A GeV. It is also supported by other experi- 
ments at RHIC energies (see references in, e.g., |j|). How- 
ever, due to the lack of data and the rather large experi- 
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mental errors, a firm conclusion is still awaiting. It is also 
urgent to perform microscopic transport-model calcula- 
tions to explore possible deviation from the ideal mr- 
scaling especially in order to understand 'non-Gaussian" 
effects as well as the HBT puzzles. 

In this paper, we investigate the my dependence of the 
HBT parameters with three kinds of identical-particle 
correlations: ir~ — ir~ , K + — K + . and A — A. So far the 
most majority of HBT investigations are on the tt — tt 
pairs due to their large yield. In order to obtain a cleaner 
signal with sensitivity to earlier stages of HICs, kaon- 
kaon correlations have been explored at several energies 
0, HH HI, HE HI| • Besides mesons, in order to map out a 
widespread m^-dependence, the A — A baryonic correla- 
tion is the next convenient choice. It is an identical non- 
charged-particle correlation so that a similarly Gaussian 
shape as for pions can be expected if the nuclear potential 
is not considered in FSI. 

To explore the m-r-scaling, we employ the UrQMD 
model. In UrQMD, the hadrons are represented by Gaus- 
sian wave packets in phase space. After the Wood-Saxon 
initialization, the phase space of hadrons is propagated 
according to Hamilton's equation of motion [3, [2^ |. 
fj = Jj^- and Pi = — f^-- Here r.; and are the co- 
ordinate and the momentum of hadron i. The Hamilto- 
nian H consists of the kinetic energy T and the effective 
interaction potential energy U, H = T + U. The two- 
body Coulomb potential is considered for formed charged 
particles. Recently, a soft equation of state with mo- 
mentum dependence (SM-EoS) for formed hadrons and 
a density dependent Skyrme-like term for "pre-formed" 
hadrons from string fragmentation have been supplied 
into the UrQMD transport model, please see details in 
(l3l I27I ]. For observables such as the nuclear stopping, 
the elliptic flow and the HBT parameters of pions visible 
improvement towards the data has been observed if these 
potentials are included. Especially, the HBT "t-puzzlc" 
for pions can be consistently solved. In this paper, we 
continue this topic and further calculate the HBT pa- 
rameters of kaon and A sources. 

Three cases of experimental data at mid-rapidity are 
compared with our calculations and the phase space cuts 
are same as those listed in Q, the three systems un- 
der investigation are: (1) central Pb+Pb collisions at 
the SPS beam energy E lab = 20 A GeV (dubbed "E20"), 
(2) central Pb+Pb collisions at the SPS beam energy 
E lab = 158A GeV ("E158"), (3) central Au+Au colli- 
sions at the top RHIC energy ^/JWn = 200 GeV ("s200"). 
For each case about 25 thousand central events are cal- 
culated. If not stated otherwise, the UrQMD transport 
program stops at t cut = 200fm/c, the residual long-lived 
unstable resonances are forced not to decay after this final 
cut time. All particles with their phase space coordinates 
at their respective freeze-out time tf (last collisions) are 
put into the analyzing program to be discussed below. 
When studying the time-evolution of the particle corre- 
lation (shown in Fig. [5]), we produce output at t cut — 5, 
10, 15, 20, 25, 50, and 75fm/c, separately, while other 



parameters are not altered. 

The "correlation after-burner" (CRAB v3.0/3) program 
[H, H^] is then adopted for analyzing the interactions of 
two particles after freeze-out with quantum statistics and 
final state modifications so that one can construct the 
HBT correlator to be compared with experimental data. 
In this work, the strong interaction of pions and kaons 
is not considered in FSI because of its negligible effect 
[301 • However, we also noticed that the effect of strong 
FSI influences the final HBT results of neutral source 
[23l ]. For As, the nuclear modification is not considered in 
this work for simplicity. We will, however, briefly discuss 
this issue in Fig. [3] Concerning the contribution of the 
final state Coulomb interaction to the correlator, it is 
checked for the charged pion-pion and kaon-kaon cases. 
One billion pairs are performed at mid-rapidity of single 
or two correlated particles as stated in [7[ in each CRAB 
analyzing run. 

In the next step, we fit the correlator as a three- 
dimensional Gaussian form under the Pratt convention 
(using ROOT [3l[ and the x-squared method), i.e., the 
longitudinally co- moving system (LCMS). When the nu- 
clear and Coulomb modifications are not considered in 
the correlator, the fitting function can be expressed in 
the standard way, 

C(q L ,qo,qs) = l + Xe-^-^o-RWs-^oLioqL (1 ) 

In Eq. (|T|), A is normally referred to as an incoherence 
factor. It might be also affected by many other factors, 
such as the contaminations, long-lived resonances, or the 
details of the Coulomb modification in FSI. Thus, we re- 
gard it as a free parameter. Rl, Ro, and Rs are the Pratt 
radii in longitudinal, outward, and sideward directions, 
while the cross-term Rol plays a role at large rapidities. 
qt is the pair relative momentum q (q = pi — P2) in the 
i direction. 

If one considers the Coulomb effect in FSI for charged- 
particle pairs, a Bowler-Sinyukov method, which has 
been used in STAR experiments [30(, can be exploited 
in the fitting process: 

C(q L ,qo,q S ) - (1 - A) + 

XK cmd (q mv )(l + e -«x^-fl 2 4-fl 2 S «|-2fl^ g o^ )( 2) 

where the K cou i is the Coulomb correction factor and 
depends only on qi nv as same as experiments. The qi nv = 
•\/q 2 — (AE) 2 is the invariant relative momentum, where 
AE = Ei — E% is the energy difference of two particles. 

Fig. [T] shows the one-dimensional invariant correla- 
tion function for pion (left plot), A (left plot), and 
kaon (right plot) sources for central Au+Au reactions 
at y/SNN = 200 GeV. A transverse momentum kx cut is 
not applied. The pairs are binned with Aqi nv — 5MeV/c 
for pions and kaons while Aqi nv — lOMeV/c for As due 
to the small yield of As in each event. The upper limit 
is qinv — 120MeV/c. The pairs within q < 5MeV/c are 
not used for fitting due to the large errors from split and 
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merged tracks in the experiments. The left plot shows 
that the Gaussian parameterization is still suitable for 
pion-pion and A-A correlators if strong and Coulomb 
modifications are not considered in FSI. It was found that 
the Coulomb correction factor influences the correlator 
of pions at qi nv < 30MeV/c [3, [S2] which is supported 
by the present calculation as indicated by the flat K^ oul 
(thin dotted line). It is also found that the K^ oul factor 
in Eq. ((2]) reduces the final HBT radii of pion source only 
slightly at < lOOMeV/c, while the reduction of A is 
about 0.1. A similar effect of Coulomb modification on 
the HBT radii and the A factor was also observed in pre- 
vious calculations [HI, [HI • The one-dimensional fitting 
results (Rinv) for pions and As show that the pion's ho- 
mogeneity length is about twice larger than for A's, which 
implies that more As are emitted from the early stage. In 
addition, the larger A factor indicates that As are less af- 
fected by the decay of resonances at late stage, although 
some amount of As are produced by the decay of S* res- 
onances. Fig. Q] (right) shows the Ri nv of the K + source 
which lies between that of the pion and A source, also 
the A value of the kaon source is largest, which supports 
that kaons are produced earliest among the investigated 
particles |25j. However, for kaons, the Coulomb modifi- 
cation alters the correlator up to the large value of qi nv , 
~ 60MeV/c. Although a strong reduction of the Ri nv is 
not seen, A is reduced by ~ 0.3. The Coulomb corrected 
A values of pions and kaons match the data much better 
f22l . [30l |35[ | . Due to the fact that the Coulomb modifi- 
cation to the kaon-kaon correlation is larger than for the 
pion-pion correlation, we adopt the Bowler-Sinyukov fit- 
ting method expressed in Eq. {2]| after constructing the 
3-dimensional correlator of kaon source with Coulomb 
modifications. 

The HBT par ameters can be as a function of mr, 
where Tot = y/k^ + m 2 where m is the pion, kaon, 
or A mass and = (pit + P2t)/2 is the transversal 
component of the average momentum k of two particles, 
k = (pi + P2V2. Fig. [5] depicts the calculated mr- 
dependence of the HBT radii Rl (top plots), Ro (middle 
plots), and Rs (bottom plots) of ir~ (dashed lines) and 
K + (dash-dotted lines) sources in central Pb + Pb colli- 
sions at Ei ab = 20A GeV (left plots) and at Ei ab = 158A 
GeV (right plots). The cascade results are compared 
with data (solid stars for pion data at both beam en- 
ergies [36j], open symbols for kaon data at 158 GeV from 
the NA44 and the NA49 collaboration [U Hi]). One 
should bear in mind that the kaon data are obtained un- 
der slightly different physical cuts from the pion case. For 
instance, in Ref . [24[ , the 5% most central interactions are 
selected for the central kaon data sample while the 7.2% 
most central interactions are selected for pions (36|. At 
Ei a b — 20A GeV, we observe m^-scaling for pions and 
kaons in Rl and Rs, only the calculated Ro of the pion 
source is slightly larger than the data. It is also seen 
that the m^-dependence of Ro of kaon source is weak. 
At Ei ab = 158A GeV, the m^-scaling of Rs and Rl still 
hold with similar problems for Ro- However, one should 
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FIG. 1: One-dimensional invariant correlation function for 
pion, A, and kaon sources for "s200" case. The potentials for 
both "pre-formed" and formed particles are considered in cal- 
culations. In the left plot: The squares and circles represent 
the calculation results for pion and A source without modifi- 
cations on FSI, separately. The dashed and dash-dotted lines 
are the fitting results to them. In the right plot: The squares 
and circles represent the calculation results for kaon source 
with and without Coulomb modification, separately. The 
dashed and dash-dotted lines are the fitting results to them. 
The Coulomb correction factors of 7r~ — 7r~ and K + — K + 
pairs (KZoui an d K^ oul ) are also shown in both plots by thin 
dotted lines. 



note that the few currently available kaon data can not 
sufficiently evaluate the performance of the scaling. 

Let us move on to see the m T -dependence of HBT 
radii of pions (dashed lines), kaons (dash-dotted lines), 
and As (dash-dot-dotted lines) at the top RHIC energy 
\J s nn = 200 GeV, which is shown in Fig. [3] The 
pion data (solid stars) are from [3(J [35|, and the pre- 
liminary kaon data (open stars) are from [22||. The 
left plots show the calculations with the SM-EoS for 
formed baryons, while the calculations with potentials 
for both "pre-formed" and formed particles are shown in 
the right plots. In each plot, the m^-scaling function 
Rl = Ro = Rs = 3/y / my is also shown by solid line. 
In Fig. [3] (left), one can not observe mT-scaling in all 
HBT directions if the potentials are considered only for 
formed baryons. Although the calculated Rs of kaons 
stay on the scaling line, those of As are about 1.2 fm 
above the scaling line. Further, Rl values of kaons and 
As are larger than the scaling results. In the Ro direction 
the same is true for the HBT radii of pions, kaons and As. 
Let us now consider potentials for "pre-formed" hadrons 
as discussed in [l3J|. Fig. [3] (right) shows that the trans- 
verse radii Ro and Rs of the pion source nicely follow 
the scaling line. However, a steeper mr-dependence of 
Rl is seen which might be due to the absence of a proper 
collision term for the "pre-formed" hadrons. Secondly, it 
is exciting to see that the "pre-formed" hadron potential 
leads to a smaller Rl (Rq) of kaon pairs as well so that 
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FIG. 2: mr-dependence of the HBT radii of -k~ (dashed lines) 
and K + (dash-dotted lines) sources in central Pb + Pb col- 
lisions for "E20" (left plots) and "E158" (right plots) cases 
(with the cascade mode in UrQMD calculations). The ex- 
perimental NA49 and NA44 data of pions (for "E20" and 
"E158") and kaons (for "E158") are also shown by solid and 
open symbols, separately [2ll. [24 l36l] . 

the HBT radii of kaons also follow the scaling line quite 
well. Thirdly, the results for As also approach towards 
the scaling line especially at large kr where the effect of 
resonance decay is minor because of an early emission. 
At small kx, the decay of long-lived resonance £(1385) 
into As enlarges the HBT radii. Meanwhile, the absence 
of the nuclear modification on FSI might also play a role, 
which deserves a further investigation. 

Based on discussions in Figs. Q] - [3j we conclude that 
interactions at the late stage (including the modifica- 
tions in FSI) indeed influences the HBT parameters. 
The idea of "pre- formed" hadron interactions does not 
only solve the HBT "time (i?o /-Rs)-puzzle" , but also 
improves the mT-scaling considerably. To elaborate fur- 
ther on this point, Fig. 2] shows the time evolution of 
the ratio Ro/Rs and the incoherence factor A of pion 
pairs at 250 < kx < 350MeV/c for central Au+Au reac- 
tions at y/sjyN = 200 GeV. Calculations with and with- 
out "pre-formed" hadron potential are compared. The 



FIG. 3: m/r-dependence of the HBT radii of n~ (dashed 
lines), K + (dash-dotted lines), and A (dash-dot-dotted lines) 
sources in central Au + Au collisions for "s200" case with 
potentials only for formed baryons (left plots) and with po- 
tentials for both formed baryons and "pre-formed" particles 
(right plots). The experimental RHIC data of pions and kaons 
are shown by solid and open stars separately 22, 30, 35]. The 
function Ri = 3/^/mr is also shown by solid line in each plot. 



corresponding experimental data of the Ro/Rs an d A 
are indicated by a square and a circle at the right end 
of the plot. Without "pre-formed" particle interactions, 
the Ro I Rs ratio increases rapidly up to t ~ 25fm/c from 
~ 1.0 to ~ 1.5. In contrast, considering the "pre-formed" 
hadron potential, the ratio Ro/Rs is nearly time inde- 
pendent (stays at its start value of Ro/Rs ~ 1 ± 0.05). 
The effects of FSI continue to play visible roles (although 
relatively weakly after 25fm/c) on the A factor, which is 
also consistent with the visible suppressing effect of the 
Coulomb modification of FSI on the A value from the 
analysis of Fig. [TJ We also find that the "pre-formed" 
hadron potential reduces the A further to approach the 
data. 

For central collisions, the HBT radii can be (approxi- 
mately) analytically obtained under the assumptions of 
thermalization and Gaussian-source shape and be ex- 
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pressed as [iH.l37l|. 
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(3b) 
(3c) 



Here the space-time coordinates x, y, z, and t are relative 



distances to their "effective source centers" (x^ 



•)): 



x n _ x fj. _ x n ^nd ^ anc j ^ T are components of the 
v elocity of particle pair (3 {(3 = k/k°, k° = (E 1 +E2)/2 ta 
\J m 2 + k 2 is the average energy of two particles. Usually, 
the on-shell approximation is used). Eq. (JHJd) can also be 
expanded as 



Rn 



(4) 



In central collisions, because of the longitudinal reflec- 
tion symmetry, (x 2 ) ~ (y 2 )- Therefore, by comparing 
Eq. ([3fc) with Eq. Q, it is clear that the difference of 
Ro and Rs mainly comes from the relative strength of 
the time- related term (/?§-P) and the x-t correlation term 
—2(PTxt) . Thus, one obtains a clear interpretation of the 
large reduction of the Ro/Rs ratio (even Ro/Rs < 1 
may happen), if "pre- formed" particle interactions are 
taken into account by relating it to the stronger phase- 
space correlation induced by the potentials. 

Fig. [5] shows the mr-dependence of the ratio Ro/Rs 
of 7r~, K + , and A sources at transverse momenta 
250 < kr < 350MeV/c for central Au+Au reactions 
at y/SNN = 200 GeV. The calculations with potentials 
only for formed baryons (lines with open symbols) and 
with potentials for both formed and "pre- formed" par- 
ticles (lines with solid symbols) are compared with pion 
data [3(| [HI and preliminary kaon data '22]. Consider- 
ing the uncertainties from the non-Gaussian effect [l6| 
as well as the variant treatments of Coulomb correction 
on FSI (see [3(3, H2J ) , a 20% deviation margin is marked 
in light gray. It is seen that, although the effect of the 
formed hadron potentials on the ratio becomes weaker 
with the increase of particle mass, the Ro/Rs values of 
pions and kaons at large kx are above the upper limit of 
the shadowed area. The inclusion of "pre- formed" par- 
ticle interactions cures these deviations and allows for 
a consistent understanding of the data. Here it should 
be addressed that due to the x-t correlation shown in 
Eq. ([4]), the Ro/Rs ~ * 1 does not mean that the mean 
duration time should approach unity, a finite non-zero 
mean proper duration time has been discovered by re- 
cent experiments |18| via three-dimensional source imag- 
ing methods which definitely deserve further transport- 
model analysis. 

To summarize, the mj dependence of the HBT pa- 
rameters for pion, kaon, and A sources are investigated 
at three SPS and top RHIC energies with a microscopic 
transport model (UrQMD). In the cascade mode (i.e., 
without potential interactions), the HBT radii Rl and 
Rs fall roughly with the mr-scaling only at the lowest 
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FIG. 4: Time evolution of the ratio Ro/Rs and the A factor 
of pion source at 250 < kr < 350MeV/c, which are calcu- 
lated with potentials for formed baryons and with potentials 
for both formed and "pre-formed" particles, are shown by 
different lines. The Ro / Rs and A experimental data are rep- 
resented by a square and and a circle at the right end of the 
plot. 



SPS energy Ei ab = 20A GeV, while the R at all en- 
ergies and the Rl at energies Ei a f, — 158 A GeV and 
\J snn = 200 GeV are apart from the mj- scaling line. 
Even the inclusion of potentials for formed hadrons does 
not alter this picture. However, with a density dependent 
potential for "pre-formed" particles from string fragmen- 
tation, the HBT radii of pion and kaon sources, even 
those of the A source, are seen to follow the mT-scaling 
line R = 3/^rriT fairly well. The strong and Coulomb 
interactions at late stage are observed to influence the 
HBT radii and the A factor. The HBT "t-puzzle" in 
many calculations exhibited by the larger Ro/Rs ratios 
of pion and kaon sources than data can only be solved 
by the "pre-formed" hadron potential which dominates 
at the early stage of HICs. 
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